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We describe the development of the model for interstellar 7 -ray emission that is the standard 
adopted by the LAT team and is publicly available. The model is based on a linear combination of 
templates for interstellar gas column density and for the inverse Compton emission. The spectral 
energy distributions of the 7 -ray emission associated with each template are determined from a fit 
to 4 years of Eermi-LAT data in 14 independent energy bins from 50 MeV to 50 GeV. We fit those 
distributions with a realistic model for the emission processes to extrapolate to higher energies. We 
also include large-scale structures like Loop I and the Eermi bubbles following an iterative procedure 
that re-injects filtered LAT counts residual maps into the model. We confirm that the cosmic-ray 
proton density varies with the distance from the Galactic center and find a continuous softening of 
the proton spectrum with this distance. We observe that the Eermi bubbles have a shape similar to 
a catenary at their bases. 


I. PRINCIPLE 


This paper describes the model of interstellar emis¬ 
sion recommended for point source analyses of the 
Fermz-LAT Pass 7 reprocessed data (P7REP) [T]. The 
high-energy interstellar 7 -ray emission is produced by 
the interaction of energetic cosmic rays (CRs) with 
interstellar nucleons and photons. The decay of sec¬ 
ondary particles produced in hadron collisions, the in¬ 
verse Compton scattering of the interstellar radiation 
field (ISRF) by electrons and their bremsstrahlung ra¬ 
diation emission in the interstellar medium (ISM) are 
the main contributors to the Galactic diffuse emission. 
The interstellar emission model is based on a template 
method: we assumed that the diffuse 7 -ray intensity 
at any energy can be modeled as a linear combination 
of maps or templates of hydrogen column-density {Ih) 
and a predicted inverse Compton (ICp) intensity map 
(IiCp) calculated by GALPROP [^. The intensity of 
each template in each energy bin is determined from 
a fit to the P7REP Clean class events from the first 4 
years of the mission binned into 14 equal logarithmic 
intervals from 50 MeV to 50 GeV. 

In addition to the interstellar emission, the LAT de¬ 
tects 7 rays from other sources that need to be taken 
into account in the analysis. We do this by adding 
dedicated components to account for a residual inten¬ 
sity of the Earth limb {Iiimh)i for point and extended 
7 -ray sources {Iext)i and for the emission from the Sun 
and the Moon {Isun-moon)- Einally we add an uniform 
intensity template {liso) to account for unresolved 7 - 
ray sources and CR contamination in the data. Eor 
a given counts map pixel and energy band we calcu¬ 
lated the predicted number of counts {Np^ed) detected 
by Fermi-LAT as: 
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Npred(E)= qi{E)lH, + Nic(E)IiCp{E) 

i=H templates 

TAtiso{E)Iiso T-^i/oopi(-F)/x,oop/ T E Atpatchi {E)Ipatchi 

i=patch 

T (-F')hi7TT,b + E Npt^(E)S(i) 

i=point src 

+ E Eexti{E)Iexti T I sun-moon {E) 

i=extend src 

( 1 ) 


where E is the energy. We use the notation I to 
denote predicted counts maps resulting from the 
convolution by the LAT PSE of the product of an 
intensity map I and the instrument exposure and 
pixel solid angle. In Equation each template of 
hydrogen column density Ih is multiplied by its 
associated hydrogen 7 -ray emissivity q. The factors 
A’/c, Niimi), and Niso represent the renormalization 
factors associated to lie pi ^Hmbi Q^nd liso- Fun-moon 
is kept fixed in the analysis. Equation also incor¬ 
porates coefficients associated to extended sources 
{Next) and to point sources {Npt) represented by 
the Dirac 5 function. Foopi and N^oopi account 
for local IC emission from Loop l. Eor unmodeled 
excesses, we also introduced in Equation patches of 
uniform intensity {Ipatch) with normalization factors 
Npatch- The procedure behind the construction of 
our interstellar emission model was to find templates 
for the gas column density and ICp, fit Equation 
to Fermi-LAIL counts maps with g, A/c, Nisoi 
NloopIi Npatch 1 Niimbi Npt and Next left free to vary 
in each energy bin and to extrapolate the coefficients 
related to the hydrogen templates and ICp outside 
the energy range of the fit. 
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II. TEMPLATE DESCRIPTION 

About 99% of the ISM mass is gas and about 70% 
of this mass is hydrogen. The hydrogen gas exists in 
the form of neutral atoms in cold and warm phases, in 
the form of neutral molecules (H 2 ), and in an ionized 
state. Helium and heavier-elements are assumed to be 
uniformly mixed with the hydrogen. 

H I is traced by its radio 21 -cm line radiation; we 
derived its column density AH l from the 21 -cm line 
radiation temperatures under the assumption of a uni¬ 
form spin (excitation) temperature (Ts). The 21-cm 
all-sky Leiden-Argentine-Bonn (LAB) composite sur¬ 
vey of Galactic H l [ 2 ] is used to determine the all-sky 
distribution of AH l. We derived AH l from the ob¬ 
served brightness temperature using a Ts of 140 K 
which provided the best fit to the Fermi-LAT data in 
regions with 90° < I < 270° and \b\ < 70°. Because 
the CR flux varies with Galactocentric distance and 
Equation is valid only if the CR flux is uniform in 
each template, we partitioned the Galaxy into Galac¬ 
tocentric annuli and assign to each annulus the corre¬ 
sponding H I column density. 

The molecular hydrogen which does not have a per¬ 
manent dipole moment, generally can not be observed 
directly in its dominantly cold phase. The observa¬ 
tion of molecular gas relies on other molecules and 
especially on the 2 . 6 -mm J= 1^0 line of carbon 12 
monoxide (CO). It is common to assume that the 
H 2 column density is proportional to the velocity- 
integrated CO brightness temperature IT(CO). The 
molecular hydrogen-to-CO conversion factor is ex¬ 
pressed as Xco=A'{}l 2 )/W{CO). We obtained the 
IT(CO) spatial distribution from the Center for As¬ 
trophysics composite survey [3]. We derived Galacto¬ 
centric annuli from radial CO velocities in a similar 
way as for H l. 

Unfortunately CO is not a perfect tracer of H 2 . More¬ 
over AH I derived under the hypothesis of a uniform 
Ts=l4:0 K is likely to be biased in lines of sight cross¬ 
ing regions of different Ts . Those approximations lead 
to large underestimates of the quantities of gas called 
dark neutral medium (DNM) in our Galaxy [4]-[6]. 
Since dust is well mixed with gas, we accounted for 
this gas by including in our model a template related 
to the total dust column density. We derived a DNM 
template as the residual map obtained after subtract¬ 
ing from the dust optical depth map of the [7] parts 
linearly correlated with the AH l and lU(CO) annuli. 
Subtracting the correlated parts from the dust opti¬ 
cal depth map revealed coherent structures across the 
sky both in the positive and in the negative residuals. 
The negative residuals are likely related to regions in 
which an average Ts of 140 K is too low, and thus 
AH I is overestimated. In this paper we call this map 
the “AH I correction map”. The positive residuals 
reveal gas in addition to that traced by AH l and 
lU(CO). Here we associated this excess map to the 


DNM distribution even if it also includes regions in 
which an average Ts of 140 K is too high and poten¬ 
tially incorporates also ionized hydrogen that might 
be mixed with dust. 

Due to their proximity, CR protons and electrons in¬ 
teracting with the Earth limb make the Earth by far 
the brightest 7 -ray source in the sky [ 8 ]. The Fermi 
standard observational strategy is such that the Earth 
is not directly in the field of view of the LAT. However 
a large number of limb photons entering the LAT at 
large zenith angles are still detected. We constructed 
a simple template based on the subtraction of the 
counts map derived with a zenith angle cut at 100 ° 
to a counts map restricted to angles above 80°. 

In fitting the model for interstellar diffuse emission we 
included templates for 21 extended sources and 2179 
point sources at positions listed in a first iteration of 
the third Fermz-LAT catalog (3EGL) [9] derived with 
a preliminary iteration of the Galactic diffuse emis¬ 
sion model. We also incorporated in the fit the 7 -ray 
emission from the Sun and the Moon. Their intensities 
were not allowed to vary during the 7 -ray template fit 
procedure. 

While the different gas column-density maps offer 
templates for photons originating mainly from tt^- 
decay and Bremsstrahlung emission, there is no direct 
observational template for the IC emission. Instead it 
must be calculated. Eor that we used the prediction 
from the GALPROP code with GALDEE identifica¬ 
tion ‘^A^6^30^ 150^2, a representative diffusive reac¬ 
celeration model described in uni. 

Excesses originating for example from Loop l [iniiii] 
or the Fermi Bubbles [Tans] are observed when we 
compared a preliminary template model derived only 
from the templates mentioned above to the Fermi- 
LAT observations. There is no accurate a priori tem¬ 
plate for the 7 -ray emission of those large structures. 
Not including them in Equation as well as other 
structures that we observed in the residuals at lower 
latitudes, will strongly bias the fit. To reduce this 
bias we roughly modeled the strongest emitting re¬ 
gion of Loop I with a selected region of the 408 MHz 
radio continuum intensity from the survey of [14] . To 
account for excesses not correlated with radio tem¬ 
plates we introduced ad hoc patches in Equation 
The patches are regions of spatially uniform inten¬ 
sity whose shapes encompass regions with an excess 
of photons of at least about 20 % compared to the pre¬ 
liminary model. We added 4 patches including a large 
rounded shape filling Loop l and three smaller patches 
closer to the Galactic plane. Additionally, we have 
created two patches for the Fermi bubbles. We also 
made a disk-shaped patch around the Gygnus region. 
We used the patches to derive the 7 -ray emissivities 
of the hydrogen that we deduced from the fit of Equa¬ 
tion to the Fermz-LAT data. But the patches do 
not provide an accurate enough description of the in¬ 
terstellar emission to be added to the final interstellar 
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model. Instead we have incorporated a filtered 7 -ray 
residuals map obtained with the patches removed as 
described below. 


III. GAMMA-RAY FIT AND EMISSIVITY 
INTERPRETATION 

We did not fit Equation to the Fermi-LAT all¬ 
sky at once, we applied latitude and longitude cuts to 
define subsets corresponding to regions where some 
templates are prominent. We fitted the local atomic 
hydrogen in the whole longitude range but away from 
the Galactic plane at latitudes \b\ > 10°. We lowered 
this latitude cut for local CO annulus to 4° and to 3° 
for the AH l correction maps and DNM template ob¬ 
tained from the negative and positive dust residual. 
For the annuli with Galactocentric radii larger than 
10 kpc (H I and CO outer Galaxy annuli) we applied 
only a longitude cut corresponding to 90° < I < 270°. 
For those 4 independent fits we left all the template 
normalization coefficients of Equation free to vary 
in each of the 14 energy bins except for the Sun and 
the Moon templates. 

The inner Galactic region is particularly difficult to 
model. The gas column densities are affected by opti¬ 
cal depth correction, self-absorption of H l and limited 
kinematic distance resolution at low longitudes. Addi¬ 
tionally 7 -ray point and extended sources are numer¬ 
ous, a DNM template is lacking, and the IC^ morphol¬ 
ogy is uncertain. Possibly for one of those reasons, or 
because of an excess of CR or an incorrect modeling of 
a foreground emission, we observed a broad unmod¬ 
eled emission (referred to as “extra emission” in the 
rest of the text) in the direction of the inner Galaxy 
with a maximum in the first Galactic quadrant at a 
longitude of ^30°. Up to this stage this emission was 
approximately accounted for by patches of uniform 
intensity. At this point we removed the patches in 
Equation!^ and we modeled the extra emission with a 
two-step iterative procedure: 

In the first step we fitted Equation[^to the Fermi-LAT 
observations excluding 5° < / < 90° and | 6 | < 20°. 
We obtained a residual map with some emission not 
accounted for, we selected the positive residuals, we 
smoothed them with a 2 -dimensional Gaussian sym¬ 
metric kernel of 3° FWHM and re-injected them in 
Equation as a template for another iteration. We 
iterated three times up to the point where the pos¬ 
itive residual intensities approximately equalled the 
intensities of negative ones. We obtained a first set 
of 7 -ray emissivities in the inner H l and CO annuli. 
Due to the difficulty of modeling the extra emission 
we reduced the number of free templates and contin¬ 
ued the procedure with a single hydrogen template for 
each annulus: AH=AH l-l-2Xc'oW(CO). We deduced 
the Xco conversion factor as half the ratio between 
the emissivity associated to the H l template and the 


one associated to the CO one. We repeated the fitting 
procedure with the single AH for inner templates it¬ 
eratively adding the positive residuals of the fit to the 
extra emission template. In this way we obtained a 
first version of the inner template emissivity and a 
template for the extra emission and large scale struc¬ 
tures. 

For the second step we fitted the whole Galaxy in¬ 
cluding the first quadrant with all the templates pa¬ 
rameters free to vary. We added to Equation the 
template for the extra emission obtained in step one. 
We repeated the iteration and derived the final emis¬ 
sivities per hydrogen atom in the inner annuli. We 
also obtained a template corresponding to the extra 
emission and large scale structures. 

Figure shows the differential 7 -ray emissivity 
per hydrogen atom = qjAE, where AE is 
the energy bin width, for the 9 Galactocentric an¬ 
nuli and the central molecular zone (CMZ) region 
scaled to emissivity per hydrogen atom assuming 
Ac'o=0-5xl0^^ cm^ (K km s“^)“^. To derive an 
interstellar diffuse emission model at energies up to 
600 GeV we fitted the emissivities with a 7 -ray 
production model of bremsstrahlung emission and 
hadronic decay and used this model for the extrapola¬ 
tion. We fitted the differential emissivity of each annu¬ 
lus between 200 MeV and 30 GeV with a parametrized 
proton flux and a 7 -ray production cross-section based 
on m- We adopted a spectral model for CR protons 
of the form: where j3 = v/c, R is the rigid¬ 

ity of the proton and (A, Pi,P 2 ) are free parameters 
[16] and folded this proton functional with the 7 -ray 
production cross-section. Permz-LAT detects 7 rays 
resulting not only from proton-proton collisions but 
also from the interaction of heavier CR or ISM nuclei. 
We used the results of m to scale the proton-proton 
cross-section to the nucleus-nucleus cross-section tak¬ 
ing into account the ISM and CR composition. The 
7 rays detected by Eermi at energies relevant for this 
work are also produced by bremsstrahlung of electrons 
and positrons. We accounted for this contribution us¬ 
ing an electron spectral form with free parameters and 
the cross section of [18] . We fitted the bremsstrahlung 
emission together with the hadron decay component 
first to the local emissivities. We derived a proton 
functional parameter Pi and a bremsstrahlung contri¬ 
bution that we assumed constant for the other annuli. 
We then fitted the other annuli with only two parame¬ 
ters: the proton spectrum normalization A and proton 
spectral index P 2 . In Figure the emissivity resulting 
from this fit is represented by a dotted line. 

To check the validity of the 7 -ray template fitting 
procedure we studied the coherence of the resulting 
proton spectral parameters between the different an¬ 
nuli. In Figure we plotted the proton spectral in¬ 
dex P 2 versus the Galactocentric distance of the an¬ 
nulus. We observe a continuous softening of the pro¬ 
ton spectra with the distance from the center of the 
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FIG. 1: (a)-(j): The 7 -ray emissivity per hydrogen atom 
in H I and H 2 phases for the CMZ and the 9 Galacto- 
centric annuli. In order to extrapolate to higher energies 
we fitted a model based on proton density and production 
cross-section (dashed line). We applied the same proce¬ 
dure for the DNM (k) and the NU I correction (1) maps. 
We did not display in the graph points with emissivities 
lower than MeV^s~^sr~^MeV~^ nor points corre¬ 

sponding to the lowest energy bin that were also not used 
in the analysis. 



FIG. 2 : Proton spectrum spectral index P 2 versus the dis¬ 
tance from the Galaxy center. The spectral indexes are 
extracted from the fit of the 7 -ray emissivities for vari¬ 
ous hydrogen Galactocentric annuli. We did not include 
any systematic uncertainty in the spectral index error bar. 
The Galactocentric distance error bar represents the radial 
width of the annuli. For comparison we draw in solid line 
the spectral index for proton energies above 100 GeV ex¬ 
tracted from the GALPROP model ‘^T^G^SO^ 150^2 pT] . 


Galaxy. The spectral indices of protons extracted 
from the emissivities measured in the first annulus 
and the CMZ correspond to a region extending ±10° 
from the Galactic center where confusion with other 
templates or point sources is possible. It might also 
be contaminated by the soft extra emission not totally 
suppressed by the iterative fit procedure. 

Figure shows the radial distribution of proton 
density integrated above 10 GeV evaluated from the 



FIG. 3: Integrated proton flux above 10 GeV versus Galac¬ 
tocentric distance. For comparison we included the inte¬ 
grated proton flux from the GALPROP '^V^G^SO^ 150^2 
model (solid line). We did not include any systematic un¬ 
certainty in the integrated proton flux error bar. Those 
systematic uncertainties could be significant for example 
close to the GG where the integrated proton flux depends 
on the Xco ratio used in the analysis. 


present work. We observe a steep CR density increase 
around 3 kpc. Again, part of this increase can be due 
to a contamination by the extra emission of the emis¬ 
sivity for the annulus extending to ±30° in longitude. 
We also observe that the inferred CR proton density 
in the CMZ is about 4 times lower than the local one 
(about 8 times lower if we assume the same Xco as for 
the local annuli). [19] suggested a lower Xco to ex¬ 
plain the anomalously low 7 -ray production compared 
to the CO column density in COS-B. Again, caution 
should be used to interpret the proton density in the 
CMZ given possible confusion with point sources or 
with the extra emission, especially at low energies. 
In Figure we also show the CR proton density pre¬ 
dicted by GALPROP ^±^6^30^150^2. We note a 
reasonable agreement with the ones derived from 7 - 
ray observations, however beyond 5 kpc the predicted 
proton density gradient is steeper than the observed 
one [2QH23] . GALPROP also predicts a broader dis¬ 
tribution around 3 kpc. 

We applied the same fit method for the coefficients 
associated with the DNM and NU l correction tem¬ 
plates obtained from the positive and negative dust 
residual. In the case of the NU l correction template 
we left all the spectral parameters for electrons and 
protons free to obtain a better fit to the data. Fig¬ 
ure [^(k) and ( 1 ) show the emissivities associated with 
the DNM and the NU l correction templates together 
with the fit. To display in this graph the emissivities 
inferred from the optical depth map in the same units 
as those from the column density maps, we divided 
the measured and fitted emissivities by an arbitrary 
gas-to-dust ratio of 4 x 10^^ cm“^ mag“^. 

From the fit of Equation we obtained a normal¬ 
ization factor Nic continuously increasing from ^1 
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at 50 MeV to ^2 at 2 GeV and then continuously de¬ 
creasing back to at 50 GeV. Above 50 GeV the 
LAT 7 -ray statistics is low causing a correlation of 
the IG template with the isotropic emission in the fit. 
We decided to rely on GALPROP predictions for the 
extrapolation to high energies and did not apply any 
scaling to IjCp for energies above 50 GeV. 


IV. MODELING THE LARGE SCALE 
STRUCTURE AND THE EXTRA EMISSION 

To include in our interstellar emission model the 7 
rays produced by phenomena that lack templates like 
the large scale structures, we first created a conven¬ 
tional interstellar emission model based on gas emis- 
sivities and IG obtained as described above. We added 
the sources from a preliminary version of the 3FGL 
catalog, the predicted Sun and Moon intensities and 
used an isotropic and limb normalization derived from 
the local H l annuli fit. Figure (left column) shows 
the positive difference between the Fermi-LAT counts 
map and the counts map expected from this model 
integrated in three energy bands: 50 MeV -1 GeV, 1 - 
11 GeV, and 11-50 GeV. As expected, since we did not 
include in the model any large structures nor patches, 
we observe positive counts residuals including the re¬ 
gions around Loop l and the Fermi Bubbles. We also 
observe the extra emission broadly distributed along 
the plane at longitudes less than 50° and to a lesser 
extent at longitudes around 315°. We also observe an 
extended excess of counts toward the Galactic center 
at the base of the Fermi bubbles. Figure shows a 
closeup of the Galactic center region representing the 
difference between the Fermi-hNY counts integrated 
between 1.7 GeV and 50 GeV and those expected from 
the conventional model using only the determined gas 
emissivities and IG emission. To reduce the contrast 
due to the bright 7 -ray emission of the Galactic plane, 
we divided this difference by the conventional model 
counts (Figure 1^) and by the square root of this num¬ 
ber (Figure [^,c,d). The fluxes of some preliminary 
3FGL sources located in the Galactic ridge depend 
to some extent on the interstellar emission model, so 
a fraction of the interstellar emission can be incor¬ 
rectly assigned to point sources. To avoid any bias we 
did not subtract them from the counts map in Fig¬ 
ure and It- They were subtracted in Figure 
and We deduce from those plots that the bases 
of the Fermi bubbles have the form of a catenary in 
which the 7 -ray emission is enhanced. We observe 
hints that this enhancement is perpendicular to the 
Galactic plane and originates from the 7 -ray source 
located in the direction of the Galactic center. As 
pointed out by [12] the ROSAT all-sky survey [24] 
shows structures similar to the one of the Fermi bub¬ 
bles. Figure shows that the Fermi bubbles have a 
similar shape within 20 ° from the Galactic center in 



phsr'^ s cm 

7 16x10"^^ 


FIG. 4: Left column: Mollweide projection in Galactic co¬ 
ordinates of the Fermz-LAT counts map after subtracting 
point and extended sources, the limb and isotropic emis¬ 
sion, and a conventional interstellar model based on fitted 
gas emissivities and scaled ICp only. The residual map 
is shown for three energy bands: 50 MeV-1 GeV (top), 
1-11 GeV (middle), 11-50 GeV (bottom). We smoothed 
those three maps with a 2-dimensional symmetric Gaus¬ 
sian of 3° FWHM. Right: intensity of the modeled large 
scale emission RE Sic at energies: 204 MeV (top), 3.4 GeV 
(middle), 22 GeV (bottom). All the maps are displayed 
with a square root scaling and a pixel size of 0.25°. 


X and 7 rays. We note that the X-ray detected by 
ROSAT are strongly absorbed at absolute latitudes 
of less that 2 ° which can produce the same artifacts 
as the use of an overestimated hydrogen emissivities 
for the inner annuli in the Fermi-LAT residual maps. 
Other extended excesses are present along the Galac¬ 
tic plane including in the Cygnus region. We did 
not observe strong negative residuals except in the 
direction of the Garina arm tangent where the model 
largely over-predicts the observations. 

We chose to model those large scale and extra 7 - 
ray emissions by assuming they all originate from IG 
interactions of a population of CR electrons with the 
cosmic microwave background radiation. The physi¬ 
cal motivation behind the production of those 7 rays 
in the energy range between 50 MeV and 50 GeV is 
not relevant if the resulting model is consistent with 
the data. For simplicity we favored a unique ISRF 
and 7 -ray production process. We fitted in each pixel 
an electron spectral form so that the conventional 
model added to this IG emission reproduces the total 
counts of the pixel. In order to reduce in the mod¬ 
eled IG map the number of undetected point sources 
and small extended structures coming from the counts 
map, we transformed the spatial distribution of the 
electron spectral parameters into wavelets and filtered 
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FIG. 5: Closeup of a region within 20° of the Galactic 
center showing the Fermz-LAT counts map integrated be¬ 
tween 1.7 GeV and 50 GeV after subtracting a conven¬ 
tional model of interstellar emission made only of emis¬ 
sion correlated with the gas and the IGp. To reduce the 
contrast between the Galactic plane emission and higher 
latitudes, we display the residual in fractional units, di¬ 
viding the residuals by the model (a), and in units of 
standard deviation, dividing the residuals by the square 
root of the model (b). In (c) we additionally subtracted 
the point and extended sources from a preliminary 3FGL 
list. In (d) we show the residual map with the large- 
scale emission modeled by RESic subtracted, it con¬ 
tains features smaller than the angular scale included in 
RES_IG. The red dashed lines correspond to the cate¬ 
nary function 10.5 x {cosh{{l — 1)/10.5) — 1) (north) and 
—8.7 X (cos/i((/ +1.7)/8.7) — 1) (south) that reproduce ap¬ 
proximately the edge of the bubbles for latitude below 20°. 
We smoothed the 4 maps with a Gaussian of 1° FWHM. 


out scales smaller than 2 °. We created a spatial mask 
drawn by hand that encompasses only regions where 
the large-scale structures and extra emission is the 
largest. We call this masked and filtered emission 
RE Sic and show it in the right column of Figure 
at energies approximately corresponding to the geo¬ 
metric average of the energy interval used to display 
the counts map of the left column. We observe a good 
agreement between the large-scale structures observed 
in the counts map and the one we have modeled. In 
Figure we show the residual map integrated be¬ 
tween 1.7 GeV and 50 GeV in the direction of the 
Galactic center when RESjc is included in the model. 
The residual map is flat apart from some emission to¬ 



-20% 0 20% 40% 60% 80% 0 200 400 600 800 


FIG. 6: Gomparison between the Fermz-LAT residuals 
shown in Figure and the ROSAT all-sky survey m 
for energies 0.73-2.04 keV. The ROSAT observations are 
displayed in 10“® counts s“^ arcmin“^. The second row 
shows the same plots as in the first row together with the 
catenary curves given in Figure 


ward the Galactic center and the Fermi bubbles. This 
emission corresponds to the small scales filtered out in 
the wavelet decomposition. 

V. RESULTING MODEL OF GALACTIC 
INTERSTELLAR EMISSION 

We derived a final model for the interstellar emis¬ 
sion from the sum of the modeled differential gas 7 -ray 
emissivities the renormalized hcp^ and from 

the large-scale and extra emission RESjc (Equation 

i- 

I{E,l,b)= Y. ^^(^E)lH,(l,b) 

i=HI,H2,DNM ( 2 ) 

^Nic{E)IiCp (E, /, b) + RESiciE, /, b) 

We compared the Fermi-LAT counts map inte¬ 
grated between 360 MeV and 50 GeV (Figure top) 
to the one predicted by the interstellar emission model 
given by Equation combined with the ones origi¬ 
nating from non-Galactic interstellar origin (Figure]^ 
middle) . We derived the residual map (Figurebot¬ 
tom) by subtracting the model from the data and nor¬ 
malizing by the square root of the model to enhance 
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FIG. 7: Top: All-sky Mollweide projection for 4 years of 
Fermi-LAT 7 -ray counts in the 0.36-50 GeV energy band. 
Middle: counts prediction in the same energy range based 
on the interstellar emission model combined with modeled 
point and extended sources (including the Sun and the 
Moon), the Earth limb emission and the isotropic emis¬ 
sion. Both maps are displayed with square root scaling to 
enhance emission away from the plane. Bottom: residual 
map in units of standard deviations after smoothing with 
a Gaussian of 2° FWHM. The pixel size for the three maps 
is 0.25°. 


deviations above statistical fluctuations. The over¬ 
all agreement between observations and model is very 
good, partly because some excesses we observed were 
modeled and re-injected into the interstellar model. 

We derived from Equation a model of the inter¬ 
stellar emission available at the Fermi Science Sup¬ 
port Center (FSSC) website as a FITS file named 
glLiem-vOO.fit. We resampled all the maps to a 0.125° 
grid. The FITS file comprises 30 logarithmically- 
spaced energy bins between 50 MeV and 600 GeV. 
It gives the differential intensity of the Galactic 


diffuse emission model in photons sr“^ s“^ cm“^ 
MeV“^. This model tuned to LAT data is not 
corrected for the energy dispersion, it can then be 
used directly with LAT data. The Pass 7_V15 IRFs 
(P7REP_CLEAN_V15) are the recommended set for 
Pass 7 reprocessed data and this model. The differ¬ 
ence with the Pass 7_V10 used for this fitting mainiy 
resides in an improved Monte Cario PSF and in 
an updated fitting procedures to determine the pa¬ 
rameters for the LAT effective area representation. 
Those minor differences modify the exposure but not 
the reconstructed LAT events. The minimum ra¬ 
tio of exposure maps (V15/V10) is 0.98 at 50 MeV 
and the maximum 1.05 at 1 GeV. In order to use 
the modei with the finai IRFs we rescaied the in¬ 
tensity by the ratio of the exposure maps evaiuated 
for each of the 30 energy bins of the Gaiactic dif¬ 
fuse emission modei. The modei is then intended 
for use with the instrument response functions ver¬ 
sions P7REP_SOURCE_V15, P7REP_CLEAN_V15, 
and P7REP_ULTRACLEAN_V15. 
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